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Abstract 

We provide an up-to-date analysis of the parameter space of the minimal supergravity 
model (mSUGRA). Novel features include the new central value of the top quark 
mass, an improved calculation of the masses of the supersymmetric particles and the 
neutral Higgs bosons, constraints from b — * s£ + £~ decays, and a careful treatment of 
the most important experimental and theoretical uncertainties. In addition to the by 
now traditional plots of the allowed region in the (mo,m^ 2 ) plane, we show allowed 
regions in the planes spanned by pairs of physical sparticle or Higgs boson masses. 
Moreover, we search for the minimal allowed masses of new particles for various sets 
of constraints. We find that in many cases the direct experimental limits from collider 
and Dark Matter searches can be saturated even in this minimal model, and even after 
including the by now quite restrictive constraint on the Dark Matter relic density. 
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1. Introduction 



The minimal supergravity model (mSUGRA) PJ Ej remains the most widely studied im- 
plementation of the minimal supersymmetric extension of the Standard Model (MSSM). It 
shares the virtues of a stable gauge hierarchy (for sparticle masses not much above a TeV) 
possible Grand Unification of all gauge interactions [3J, and a plausible Dark Matter 
(DM) candidate jSJ |H| with all variants of the MSSM 2 . Moreover, it manages to describe 
phenomenologically acceptable spectra with only four parameters plus a sign: 

m , TO1/2, A), tan/3, sign/z. (1) 

Here m is the common soft supersymmetry breaking contribution to the masses of all scalars, 
mi/2 the common supersymmetry breaking gaugino mass, and A the common supersymme- 
try breaking trilinear scalar interaction (with the corresponding Yukawa coupling factored 
out); these three soft breaking parameters are taken at the scale Mx of Grand Unification, 
which we define as the scale where the properly normalized SU{2)l and U{1)y gauge cou- 
plings meet. Finally, tan/3 is the ratio of the vacuum expectation values (vev's) of the two 
Higgs doublets at the weak scale, which we identify with the geometric mean of the soft 
breaking stop masses, and /i is the supersymmetric higgs(ino) mass parameter. 

It should be admitted that the choice of parameters Q is not particularly natural from 
a theoretical point of view: why should the scalar masses and trilinear A parameters all 
be exactly the same exactly at scale Mx? From the perspective of supergravity theory, 
universality would seem to emerge more naturally at a scale closer to the (reduced) Planck 
mass, Mp ~ 2.4 ■ 10 18 GeV, if at all. However, while the possible unification of the gauge 
interactions makes a strong argument for a "grand desert" between the sparticle mass (or 
weak) and Grand Unified scales, physics at energy scales above Mx remains very speculative. 
At least as a first approximation it is therefore not unreasonable to impose our boundary 
conditions at scale M X - 

The ansatz (0) also has important virtues, in addition to its simplicity. It allows a 
quite varied phenomenology without violating any known constraints. In particular, the 
assumed flavor universality implies that supersymmetric flavor changing neutral current 
(FCNC) effects occur only radiatively, through renormalization group (RG) evolution. This 
keeps FCNC manageable, although, as we will see, flavor changing b — > s'y and b — > s£ + £~ 
decays do impose important constraints on the parameter space. A very appealing feature 
of mSUGRA is that it implements radiative breaking of the electroweak gauge symmetry [7] , 
i.e. the RG evolution naturally drives the squared mass of one of the Higgs fields to negative 
2 A good DM candidate only emerges if R— parity is conserved, which we assume throughout. 
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values, keeping all squared sfermion masses positive. This allows to determine the absolute 
value of n as function of the other parameters. 

In spite of these successes, in recent years there has been a proliferation of analyses 
extending mSUGRA. Some of these extensions [H] are based on specific Grand Unified models, 
and thus have independent motivation from theory. However, in many phenomenological 
analyses universality between sfermion masses and/or the universality of soft breaking Higgs 
and sfermion masses is relaxed jHj without specific theoretical motivation 3 . Indeed, there 
seems to be a perception that the parameter space of the model is getting "squeezed" by 
ever tightening constraints. 

Much of this perception probably comes from the by now quite accurate determination 
of the relic density of Dark Matter (DM) particles. At least in the framework of standard 
cosmology with a more or less scale invariant primordial spectrum of density perturbations, 
the analysis of large cosmological structures allows to infer the present DM density; the 
mapping of the microwave sky by WMAP plays an especially important role here [H]. This 
translates into a quite tight constraint on mSUGRA parameter space ^2] under the standard 
assumption that all DM is formed by lightest superparticles (LSPs), which were in thermal 
equilibrium after the last period of entropy production 4 . 

It should be emphasized that this tight constraint should not be cause for alarm. After 
all, the determination of the DM density, if taken at face value, constitutes a genuine signal 
of physics beyond the Standard Model (SM). Conceptually it should thus be considered 
on a par with, say, the measurement of a selectron mass. The fact that mSUGRA can 
accommodate this measurement is a further success of this model. 

Nevertheless it seems timely to re-assess the mSUGRA model, taking recent theoretical 
and experimental developments into account. Besides the WMAP (and related) cosmological 
data, these include: 

• More accurate calculations of leading two-loop corrections to the masses of neutral 
Higgs bosons JH] , which makes it somewhat easier to satisfy the stringent Higgs search 
limits from LEP for a fixed value of the top quark mass; 

• The new, somewhat lower central value of the mass of the top quark ^3|, which in 
turn goes in the direction of decreasing the predicted mass of the lightest neutral Higgs 
boson; 

3 Bucking this trend, there have also been a few recent studies where additional relations between the 
parameters in JU are imposed 

4 In standard cosmology this means the end of the inflationary epoch. 
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• Improved limits on radiative b decays and, in particular, first information on b — > s£ + £~ 
decays, which excludes scenarios where the sign of the amplitude of b — > 57 decays is 
opposite to the SM prediction [To] : 

• A growing (though not global) consensus (THj that the SM prediction for hadronic 
contributions to the anomalous dipole moment of the muon based on data from e + e~ 
colliders is more reliable, which again elevates the discrepancy between the measure- 
ment [T7J of — 2 and its SM prediction ^H] to the level of ~ 2.5 standard deviations. 

Several analyses of this kind have appeared in the last few years [T2J, whose results 
broadly agree with our's if we take the old value of the top mass, m t = 178 GeV. The effect 
of the new, reduced (central) value of mt has so far only been analyzed in refs. [T9tl20|: these 
papers have little overlap with our's. Usually the results of mSUGRA parameter scans are 
presented as allowed regions in the (mo, m^) plane. We also present similar allowed regions 
in the planes spanned by physical sparticle or Higgs boson masses; this should give a more 
direct overview of the kind of spectra that can be generated in mSUGRA. To the best of our 
knowledge, similar results have previously only been published in the (by now quite dated) 
review article [2*Tj . 

Moreover, we put special emphasis on a careful treatment of theoretical and experimental 
uncertainties. This allows us to derive conservative lower bounds on the masses of superpar- 
ticles and Higgs bosons in mSUGRA, for different sets of assumptions. We find that in many 
cases the direct experimental search limits can be saturated even if all relevant constraints 
are taken at face value. The main exceptions are the masses of the gluino and of first and 
second generation squarks, which are forced by the assumption of gaugino mass unification 
to lie at least 100 to 150 GeV above the current Tevatron limits. Imposing the DM constraint 
does not affect these lower bounds very much. 

The remainder of this article is organized as follows. In the next chapter we briefly 
describe calculational details and the constraints we impose. Sec. 3 updates ref. [22] by 
showing the allowed regions in the (m , m^) plane, as well as in the planes spanned by 
pairs of physical masses, for a few values of tan/5. Sec. 4 is devoted to a discussion of 
minimal masses of sparticles and Higgs bosons that are compatible with various sets of 
constraints. Finally, in Sec. 5 we summarize our main results and draw some conclusions. 

2. Scanning Procedures 

We use the FORTRAN program SuSpect [23J to calculate the spectrum of superparticles and 
Higgs bosons. Since these masses are defined at the weak scale while the dimensionful input 
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parameters in are denned at the scale of Grand Unification, the program has to integrate 
the set of coupled renormalization group equations connecting these two scales; SuSpect 
now uses two-loop equations for all relevant quantities (gauge and Yukawa couplings, //, 
and the soft breaking parameters). The program also computes the one-loop and dominant 
two-loop corrections to the Higgs potential, as well as the dominant one-loop corrections 
turning the running (DR) quark, lepton and sparticle masses into on-shell (pole) masses. 
The calculation of the neutral Higgs boson masses includes leading two-loop corrections ^3] • 
See ref. [2S] for further details on the calculation of the spectrum. 

Not all combinations of input parameters lead to radiative SU(2)l x U(1)y — > £/(1)qed 
symmetry breaking. This imposes a first constraint on the parameter space. We also exclude 
parameter sets where the scalar potential has deep minima breaking charge and/or color at 
the weak scale [25] . As usual in the literature [123, we do n °t ve ^° scenarios where the 
absolute minimum of the scalar potential occurs for field values intermediate between the 
weak and GUT scales j26j . since the tunneling rate into these minima is exceedingly slow. In 
the language of ref. we impose the "CCB" constraints, which exclude very large values 
of | Aq\/ JrriQ + rn\j 2 , but do not impose the "UFB" constraints. 

We next impose experimental constraints. To begin with, the strong upper limits on the 
abundance of stable charged particles [27j exclude scenarios where the lightest superparticle 
is charged. This excludes cases with m <C where f\ tends to be the LSP (especially at 
large tan/3), and some combinations with > m\/2 and sizable \Aq\, where t\ is the LSP 5 . 

We also impose the lower bounds on sparticle and Higgs masses that result from collider 
searches. We interpret the LEP limits from searches for (unstable) charged superparticles 
as requiring 

a(e + e~ -> XX; yfs = 209 GeV) < 20 fb (2) 

separately for each relevant mode (X = ti,fi,xt)- This effectively imposes lower bounds 
of 104.5 GeV on the mass of the lighter chargino xt> 101.5 GeV for the lighter scalar top 
eigenstate £i, and 98.8 GeV for the lighter scalar r eigenstate f\. For non-pathological 
situations these agree closely with the limits published by the LEP experiments fZ7\ I30j. 6 

The limits from the searches for Higgs bosons at LEP also impose important constraints 
on the parameter space. Of special importance is the limit on e + e~ — > ZH with H — > bb. 

5 These constraints can be evaded if the LSP resides in the hidden sector (e.g., if it is the gravitino |28|L 
or in extensions of mSUGRA where the LSP is an axino |29j . 

6 The only "pathological" situation that can be relevant in mSUGRA is the case of small fi — Xi mass 
splitting. However, at small tan j3 selectron searches at LEP will lead to constraints on the parameter space 
that are nearly as strong as those from t\ searches. At high tan /?, scenarios with small slepton and neutralino 
masses are excluded by the constraint. 
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In the SM it leads to the bound [31] > 114 GeV. For small and intermediate values 

of tan/? this bound applies directly to the light scalar Higgs boson in mSUGRA, but for 
tan (3 > 50 its coupling to the Z boson can be suppressed significantly. We parameterize this 
dependence as in ref. [22], except that the constant (coupling-independent) term is increased 
by 0.5 GeV in order to reflect the increase of the limit that resulted from combining the limits 
from the four LEP experiments. 

We also include constraints from quantum corrections due to superparticles. These in- 
clude the upper bound 

5p S usY < 2.2 • 10- 3 (3) 

on the supersymmetric contribution to the electroweak p— parameter [32], including 2-loop 
QCD corrections [33J. However, it turns out that this constraint is always superseded by 
either the LEP Higgs search limit or by the CCB constraint. 

A more significant constraint arises from the precise measurements of the anomalous 
magnetic moments of positively and negatively charged muons [T7j. As well known by now, 
the interpretation of this measurement hinges on whether data from semileptonic r decays 
are used for the evaluation of the SM prediction or if one only relies on data from e + e~ 
annihilation into hadronic final states. In the former case the measurement agrees quite well 
with the SM, whereas in the latter case the prediction falls ~ 2.5<r short of the measurement 
|18j . In order to reflect this uncertainty, we impose either the more conservative constraint 

-5.7-10- 10 <a AtiSU SY<4.7-10- 9 , (4) 

which describes the overlap of the 2a limits derived from the two competing SM predictions, 
or the more aggressive constraint 

1 .06 • 10- 9 < a M , susy < 4.36 ■ 10- 9 , (5) 

which is the 90% c.l. range derived using the e + e~ data only. Here gi^susy is the sparticle 
loop contribution to a M = (g^ — 2)/2. The SM prediction based on e + e~ data is nowadays 
considered to be more reliable JH] • Note that (@J) allows the supersymmetric contribution to 
vanish, or even be slightly negative, whereas © requires it to be positive. Our calculation of 
ci^ t susy is based on ref. [31] , modified to include leading logarithmic QED 2-loop corrections 
[33] which increase a Mi susY by ~ 5%. 

The constraints discussed so far are all quite robust against minor changes of the model. 
In particular, deviating from exact universality of scalar masses or, equivalently, allowing 
small flavor non-diagonal entries of the sfermion mass matrices, will not change any of these 
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bounds significantly. This is quite different for the bounds from inclusive b — > 37 decays, 
which are also widely included in analyses of the parameter space of mSUGRA and similar 
models [HI H3 EEH U2] • Including theoretical uncertainties of the SM prediction [SB] as well 
as the experimental measurement j2Z| (now statistically dominated by BELLE data), we 
require the calculated branching ratio to fall in the range 

2.65 • 10~ 4 < B{b -> si) < 4.45 • 1(T 4 . (6) 

Our calculation of this branching ratio is based on ref. j2Z], which - for heavy sparticles - 
includes the dominant QCD corrections to the x^t loop corrections, which (together with 
tH^ loops) dominate the supersymmetric contributions in all supersymmetric models where 
flavor violation is assumed to be described entirely by the Kobayashi-Maskawa (KM) matrix. 
For large tan f3 and not too heavy sparticles these contributions can be quite large; for /1 > 
they can even flip the sign of the amplitude relative to the SM prediction, leading to a second 
allowed region [22| when the modulus of this amplitude approaches its SM value. However, 
it has recently been argued ^5] that new data on b — ► s£ + £~ data strongly disfavor this 
possibility, since such a flip of sign would change the interference between penguin diagrams 
(similar to those contributing to b —>■ sj) and (new) box diagrams. We therefore impose the 
additional constraint that the amplitude for b — > 37 decays should have the sign predicted 
in the SM. 

The range (JBJ) should perhaps be extended somewhat, since the MSSM prediction has 
larger theoretical uncertainties than that in the SM. To begin with, ref. jSZ] includes NLO 
QCD corrections to the supersymmetric contribution only in the limit of heavy sparticles, 
as remarked above. Note also that the determination of the KM element Vt s , to which all 
contributions to the amplitude describing b — > 57 decays are proportional, can be affected 
by supersymmetric contributions to processes in the K sector, which in the SM lead to 
tight constraints on this quantity. More importantly, the constraint imposed by © on 
the parameters listed in will evaporate entirely if a modest amount of b — s mixing is 
allowed at the GUT scale [2E|- This would lead to one-loop gluino-mediated box diagram 
contributions to b — > sj decays [HH]- Since for strict scalar universality all contributions are 
suppressed by a factor \V ts \ ~ 0.04, even a small amount of b — s mixing would lead to new 
contributions of comparable size. The sign of this new contribution is given by the sign of 
the corresponding mixing angle, which is a free parameter in this slightly extended model. 
At the price of modest fine-tuning one could thus make any set of mSUGRA parameters 
"6 — ► 57 compatible" . Since the required squark flavor mixing would still be quite small, it 
would have negligible effects on (flavor conserving) signals at colliders, radiative corrections 
to the MSSM Higgs sector, etc. The constraint © thus has a different status from the 
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constraints discussed earlier. In Sec. 4 we will therefore present limits on sparticle and Higgs 
boson masses with or without this additional constraint. 

The last, and very restrictive, constraint that is usually imposed in analyses of constrained 
supersymmetric models is based on the determination of the density of non-baryonic Dark 
Matter (DM) from detailed analyses of the anisotropics of the cosmic microwave background 
(CMB), in particular by the WMAP experiment HTJ. At 99% c.L, they find 

0.087 < tt DM h 2 < 0.138. (7) 

Here Q measures the mass (or energy) density in units of the critical density, whereas h 
is the scaled Hubble constant. The emergence of the cosmo logical "concordance model" is 
undoubtedly a great triumph of modern cosmology. One should nevertheless be aware that 
the result ((7|) is based on several assumptions, which are reasonable but not easy to cross- 
check independently jlOj- In particular, one has to assume that simple inflationary models 
give essentially the right spectrum of primordial density perturbations. Without this, or a 
similarly restrictive, ansatz for this spectrum, the result (J7|) would evaporate. In the absence 
of a generally accepted estimate of the theoretical uncertainty from the assumed ansatz of 
the primordial density perturbations, we decided to take the 99% c.L region of the DM 
relic density, as opposed to the 90% or 95% c.l. regions used for quantities measured in the 
laboratory. 

In order to translate the constraint (JZJ) into a constraint on mSUGRA parameters, one 
has to make several additional assumptions. The lightest neutralino Xi must be essentially 
stable, which, in the context of mSUGRA with conserved R— parity, requires it to be heavier 
than the gravitino [28]. In addition, one must assume that standard cosmology (with known 
Hubble expansion rate, and no epoch of entropy production) can be extrapolated backwards 
to temperatures of at least ~ 5% of rn^o. In that case Xi was m full thermal equilibrium 
with the SM plasma, making today's Xi density independent of the "re-heat" temperature 
of the Universe Tr at the end of inflation. With these assumptions, our calculation of the 
relic density proceeds as in ref. [22]. Without these assumptions, no meaningful constraint 
on the parameters (0) results. In Sec. 4 we therefore again present results with or without 
the constraint (J7J). 

3. Results of the scans 

3.1 (mo,mi/2) parameter space 

Examples of scans of the (mo, wii/2) plane are shown in Figs. 1 through 6. Figs. 1 and 2 
represent our "base choice", m t = 172.7 GeV (the current central value JI]), A) = and 
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fx > 0. The light grey regions are excluded by theoretical constraints (in particular, by 
the requirement of correct electroweak symmetry breaking), as well as by the searches for 
sparticles, i.e. by the constraint (J2J). The dark grey areas are excluded by the requirement 
that the LSP be neutral, in particular by rrif 1 > m%o. The pink regions are excluded by 
searches for neutral Higgs bosons at LEP. The light pink regions are excluded even if we 
allow for a 3 GeV theoretical uncertainty in the calculation of m^ in the MSSM, i.e. if we 
only exclude (SM-like) Higgs bosons with calculated mass m^ < 111 GeV. In the medium 
pink region the predicted rrih lies between 111 and 114 GeV, which might be acceptable if 
unknown higher order corrections are sufficiently large and positive. As well known, the 
LEP data [31] show some (weak) evidence for the existence of an SM-like Higgs boson with 
mass near 114 GeV; the mSUGRA regions that can explain this small excess of Higgs-like 
events are shown in red. 

In the blue regions the SUSY contribution to the anomalous magnetic moment of the 
muon falls in the range (jSJ) favored by e + e~ data; recall that in this case a positive SUSY 
contribution is required at the ~ 2.5 a level. The green region is excluded by the constraint 
(jUj) on the branching ratio for radiative b — ► 57 decays. Finally, in the black regions thex? 
relic density lies in the desired range (J7J). 

We see that the Higgs search limits are very severe for small and moderate values of 
tan/3, and/or for smaller values of m t ; for example, for m t = 172.7 GeV, tan/3 = 5 and 
A = (Fig. 2) they imply m 1( / 2 > 0.6 TeV for small m , or m > 2 TeV for small m^. 
At first the region excluded by this constraint shrinks quickly with increasing tan /3, but it 
becomes almost independent of this parameter once tan/? > 20. These regions become even 
larger if m t is below its current central value. For example, Fig. 3 shows that the region 
excluded by this constraint for tan/3 = 10 and m t = 166.9 GeV, the current 95% c.l. lower 
bound on this quantity, is very similar to that for tan/3 = 5 and m t at its current central 
value. Conversely, again for tan/3 = 10, increasing m t to the previous central value of 178 
GeV |Hj (which happens to be quite close to the current 95% c.l. upper bound) reduces the 
lower bound on m 1( / 2 for small mo to about 250 GeV, and allows values of mo down to about 
750 GeV even if TO1/2 is small, see Fig. 4. All these numbers allow for a 3 GeV theoretical 
uncertainty in m^. 

The (green) regions excluded by the b — > 57 constraints shows the opposite dependence 
on tan/3, becoming larger as this parameter increases. Note, however, that for m t = 172.7 
GeV and Aq = (Figs. 1 and 2) the Higgs constraint supersedes the b — > 57 constraint for 
values of tan j3 up to tan /3 = 50, unless one assumes that unknown higher-order corrections 
to m h are large and positive. For the highest value displayed, tan/3 = 58, the b — > s-y 
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Figure 1: The mSUGRA (mi/2, mo) parameter space with all constraints imposed for Aq = 
0,/i > and tan/? = 10. The top quark mass is fixed to the new central value, m t = 172.7 
GeV. The light grey region is excluded by the requirement of correct electroweak symmetry 
breaking, or by sparticle search limits. In the dark grey region f\ would be the LSP. The light 
pink region is excluded by searches for neutral Higgs bosons at LEP, whereas the green region 
is excluded by the b —>■ sj constraint (0|). In the blue region, the SUSY contribution to g^ — 2 
falls in the range J3J) ; whereas the red regions are compatible with having an SM-like Higgs 
boson near 115 GeV. Finally, the black regions satisfy the DM constraint 

constraint excludes an additional domain close to the region where electroweak symmetry 
breaking does not take place; in this area, the charged Higgs boson is relatively light and 
the tH contribution (which is not compensated by the tx ± ones, as the top squarks are 
rather heavy) leads to a value of B{b — ► sj) that is slightly higher than the upper bound 
5 max = 4.45 • lO" 4 . 

Figs. 3 and 4 show that reducing (increasing) m t further reduces (increases) the im- 
portance of the b — > 57 constraint relative to the Higgs search constraint. Taking sizable 
(negative) values of Aq also increases the relative importance of the b — > S7 constraint, as 
shown by Figs. 5 and 6. 7 



7 The effects of A ^ have recently been studied in 
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Figure 2: The mSUGRA (mi/2, mo) parameter space with all constraints imposed for Aq = 
0,/i > and the values tan/3 = 5 (top left), 30 (top right), 50 (bottom left) and 58 (bottom 
right). The top quark mass is fixed to the new central value, m t = 172.7 GeV. Notation and 
conventions are as in Fig. 1. 

For example, for tan (3 = 30 and Aq = —2 TeV (Fig. 6) the region excluded by the b — > 57 
constraint even excludes some combinations of parameters with predicted rrih well above 115 
GeV. The predicted branching ratio for b — >• 57 decays is quite sensitive to A since the main 
supersymmetric contribution comes from stop-chargino loops, and stop mixing is affected 
quite strongly by Aq, as long as \A \ > m^. Recall, however, that this constraint can be 
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Figure 3: The mSUGRA (mi/2,^0) parameter space with all constraints imposed for A = 
0,fi > and the values tan/5 = 10 (left) and 50 (right). The top quark mass is fixed to 
mt = 166.9 GeV. Notation and conventions are as in Fig. 1. 
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Figure 4: The mSUGRA (m^jMo) parameter space with all constraints imposed for A = 
0,/i > and the values tan/3 = 10 (left) and 50 (right). The top quark mass is fixed to 
m t = 178 GeV. Notation and conventions are as in Fig. 1. 
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more easily circumvented than the other constraints discussed here. 

The excluded grey regions also increase with increasing tan /?. Increased fz — tr mixing 
as well as reduced diagonal f masses due to RG effects imply that the requirement m fl > 
m^o excludes an increasing area with m\i 2 ^> m\. Similarly, the requirement of correct 
electroweak symmetry breaking excludes an increasing area at m 2 , 3> m\, 2 . For small values 
of tan /? the increase of this excluded region is mostly due to the reduction of the top Yukawa 
coupling, which scales like 1/ sin/?; for larger tan/? the effects of the bottom Yukawa, which 
scales like 1/ cos/? ~ tan/? for tan 2 /? 3> 1, in the RGE become more important, which to 
some extend counter-act the contributions of the top Yukawa coupling. For the large value 
tan/? = 58 (Fig. 2), this region covers most of the parameter space as the bottom Yukawa 
coupling becomes very large. For even larger values, tan/? > 58, one cannot obtain correct 
electroweak symmetry breaking. 

The grey regions also depend very sensitively on the top mass, as shown by Figs. 3 and 4. 
Figs. 5 and 6 show that increasing \Aq\ reduces the region at Wq 3> mf^ 2 where one cannot 
break the electroweak symmetry; on the other hand, a significant region near the origin is 
now excluded by the false vacuum ("CCB") constraints. 

The (blue) region favored by the measurement of — 2 (if the SM prediction using data 
from e + e~ annihilation can be trusted) also expands towards larger values of and mi/2 
as tan/3 is increased. This region depends only slightly on m t ; somewhat larger values of mo 
become compatible with this constraint if m t is reduced. This is due to the reduction of \i 
caused by reducing m t . The dependence of this region on Aq is again rather mild; however, 
for Aq = —2 TeV the entire blue region is excluded by the b — > constraint. 

Finally, for tan/? = 5 (Fig. 2), the (black) regions satisfying the DM constraint (J7J) lie 
right at the border of the theoretically allowed parameter space: the stau co-annihilation 
region, where ~ m^o, lies next to the region excluded by a charged LSP, whereas for 
small tan/? the "focus point region", where \i < M\ at the weak scale, is right next to the 
region where the electroweak symmetry can no longer be broken. The same holds true for 
tan/? = 10 and m t = 178 GeV (Fig. 4), for tan/? = 30 and A = -2 TeV (Fig. 6), and for 
tan/? = 30, Aq = —1 TeV and m t = 178 GeV (Fig. 5); in this last case the small black 
region at m ^> m 1( / 2 is allowed due to almost resonant h exchange [4*3*] . 

However, for larger tan/3 and/or smaller m t there are sizable regions of parameter space 
where the DM relic density comes out too low. This happens in particular for m 2 , 3> my 2 , 
where one can have an almost pure higgsino as LSP once tan /3 > 10 (for m t = 172.7 GeV; for 
smaller m t this region appears at smaller tan/?). Comparison of Figs. 1-2 and 3 shows that 
the lowest values of m in this higgsino-LSP region decreases very quickly with decreasing 
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Figure 5: The mSUGRA {m\/2, tuq) parameter space with all constraints imposed for A = — 1 
TeV, ji > 0, tan/3 = 30. The top quark mass is fixed to m t = 172.7 (left) and 178 GeV 
(right). Notation and conventions are as in Fig. 1. 
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Figure 6: The mSUGRA (mi/2, mo) parameter space with all constraints imposed for Aq = —2 
TeV, fi > 0,tan/3 = 30. The top quark mass is fixed to m t = 172.7 (left) and 178 GeV 
(right). Notation and conventions are as in Fig. 1. 
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m t . The reason is that larger m t imply a smaller (often: more negative) soft breaking 
contribution to the squared mass of the Higgs boson that couples to the top quark, which 
in turn implies a larger value of via the conditions of electroweak symmetry breaking. 

For the same reason, a reduced m t makes it easier to find solutions with vtla — 2m^o, 
yielding the very prominent DM-allowed 11 A— pole" region visible for tan (3 = 50 in Figs. 2-4; 
recall that /i 2 contributes with positive sign to m\. This region becomes broader with reduced 
m t since lower values of also mean stronger gaugino-higgsino mixing in the neutralino 
sector, and hence a larger XiXi^ coupling. For tan/3 = 58 (Fig. 2) is everywhere 
significantly below 2m x . The very large (A,H)bb couplings nevertheless imply that for 
moderate values of m and m^, Xi annihilation through virtual A and H exchange has 
the right strength. This DM-allowed region is connected to the A— pole region at smaller 
(although still large) values of tan/3. For tan/3 = 58 and small ttiq and m^, the virtual 
A, H— exchange contributions even lead to too small a Xi relic density; however, this region 
of parameter space is excluded by the b — > 57 constraint. 

On the other hand, a reduced top mass of 172.7 GeV also implies that the "bulk" regions, 
where the DM constraint (J7|) is satisfied due to the exchange of light sleptons in the t— and 
it — channel, now lies deep in the region excluded by Higgs searches at LEP. Fig. 4 is a 
reminder that a bulk region compatible with all constraints (with the possible exception 
of the theoretically somewhat shaky b — > 37 constraint) still exists for m t = 178 GeV and 
(sufficiently) large tan j3. Recall that increasing tan (3 increases ft — tr mixing, which in turn 
increases the S— wave LSP annihilation cross section through f exchange |4^j . 

Note finally that the additional possible region where the DM constraint could be satis- 
fied, i.e. with co-annihilation of the LSP neutralino with top squarks [45,, is disfavored in 
the mSUGRA scenario that we are discussing here. 

It is interesting to note that several indications for "new physics" can be explained 
simultaneously within mSUGRA. The reduction of the central value of m t has made it a 
bit more difficult to satisfy the (aggressive) g M — 2 requirement (jSj), which prefers moderate 
values of sparticle masses unless tan (3 is quite large, in potential conflict with the LEP Higgs 
search limits. However, if we allow a 3 GeV theoretical uncertainty in the calculation of m^, 
solutions satisfying © can be found for all tan/? > 8 for m t = 172.7 GeV and A = 0; if 
finite values of Aq are considered, the lower limit on tan (3 is reduced even further. On the 
other hand, if we take the prediction of at face value, again taking m t = 172.7 GeV we 
need tan/3 > 12 (7) for vanishing (arbitrary) A ; for m t = 166.9 GeV, as in Fig. 3, these 
lower bounds increase to 20 and 10, respectively. In all these cases we can satisfy the DM 
constraint (J2J) in the fi co-annihilation region, and have a CP-even Higgs boson near 115 
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GeV, as hinted at by LEP data, while satisfying the — 2 constraint (JHJ) at the same time. 

Fig. 2 shows that these three constraints can also be satisfied simultaneously in the 
A— pole region, if tan/3 is very large. However, we did not find any points in the "focus 
point" region where the aggressive — 2 constraint can be satisfied, if we take the prediction 
for rrih at face value. In this case increasing mt, and/or introducing non- vanishing Aq, would 
allow to satisfy the Higgs search limits for smaller values of m , thereby increasing o^susy! 
however, at the same time it would increase pushing the Xi re hc density to unacceptably 
large values. Even if we only demand the calculated to exceed 111 GeV (i.e. assume 
sizable positive higher order corrections to m/J one can only reach the lower end of the range 

in this region of parameter space. 

3.2 Parameter space with physical masses 

We now present some results for physical masses. In order to keep the number of figures 
manageable, we only show results for the central value of m t = 172.7 GeV, A = 0, and two 
values of tan (3. 




200 400 600 800 LSP 200 400 600 800 m LSP 

Figure 7: The mSUGRA (m$o, m^ R ) 'parameter space with all constraints imposed for A = 0, 
/i > 0, tan/3 = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 
Notations and conventions are as in Fig. 1, except that the light grey region now also includes 
combinations of masses that are never realized in mSUGRA. 

We begin with the (m^,m e - s ) plane depicted in Fig. 7; these masses largely determine 
the phenomenology of e# pair production at e + e~ colliders. These figures actually look 
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quite similar to the corresponding results in Figs. 1-2. The reason is that over most of 
mSUGRA parameter space, Xi is dominated by the Bino component, with m^o ~ 0Ami/ 2 , 
whereas m~ fl ~ + 0.15m^ 2 implies that the mass of cr is usually quite close to m . Note, 
however, that the excluded region at m\ ^> m~ does not grow with increasing tan /3, in 

H X.~i 

contrast to the analogous region in the (mo, ^1/2) plane. The reason is that the boundary 8 
of this region is set by the search limit for (higgsino-like) charginos at LEP, which essentially 
fixes the mass of x®, which is also higgsino-like here. 

The (mg,mu L ) plane is shown in Fig. 8; since the other first and second generation 
squarks have masses quite close to m^ L this plane essentially determines the cross section for 
the production of strongly interacting sparticles at hadron colliders jlH] (with the possible 
exception of a light t%\ see below). In this case both masses depend significantly on the 
gaugino mass parameter, with m? L ~ Wq + 6mf i % and m g ~ 2.5m!/ 2 - As a result, the 
accessible part of parameter space gets squeezed, whereas the entire region < 0.8m.g is 
not accessible |17j , greatly increasing the size of the grey regions compared to the analogous 
results of Figs. 1-2. Moreover, since m g is independent of /i, the region at rri~ lL 3> m~ that 
is excluded because /z 2 comes out too small does grow with increasing tan /3. Note that our 
basic parameter scan only included values of mo up to 3 TeV. As a result, the area with 
m^ L > 3 TeV and much smaller rrig did not get probed, although some of it is theoretically 
accessible. 

We next turn to the (m^mfj parameter space depicted in Fig. 9. In this case both 
masses depend significantly on tuq, but only m^ has a strong dependence on mi/2- The 
"focus point" region with higgsino-like or mixed LSP is therefore still at the top-left of the 
accessible region. However, the inaccessible region at the left side of the figure is considerably 
larger than in Figs. 1-2, since for our choice Aq = 0, one cannot have too much above 
mj . Note also that we chose different y— scales in the two frames of Fig. 9. Increasing tan/3 
increases the r Yukawa coupling, which reduces the soft breaking masses in the f sector 
through RG effects, and increases tl — tr mixing; both effects tend to reduce m ?1 . For 
tan/3 = 10, mfj is quite close to m e - K , but for tan/3 = 50 it is significantly smaller. On the 
other hand, m^ is relatively insensitive to tan /3. 

The Higgs sector [3B] reflects the radiative symmetry breaking in mSUGRA. For small and 

moderate values of tan j3 the heavier Higgs bosons, whose masses are essentially determined 

by that of the CP-odd Higgs boson A, are among the heaviest of all new particles [12] • On 

the other hand, for large tan j3, RG effects due to the bottom Yukawa coupling greatly reduce 

rriA ■ We saw in the previous section that this leads to scenarios where ttla — 2m^o , and hence 

8 Most of this region is excluded since electroweak symmetry breaking would require [x 1 < 0; however, 
there is also a small area where fi 2 is positive, but below the LEP limit. 
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Figure 8: The mSUGRA (nig, m^ L ) parameter space with all constraints imposed for Aq = 0, 
a > 0,tan/3 = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 

to strong Xi annihilation in the early Universe due to near-resonant A boson exchange. This 
reduction of is reflected by the reduced y— scale in the right frames in Figs. 10 and 12, 
which show the (m^o, m^) and (rrih, niA) planes, respectively. In fact, Fig. 10 shows that for 
tan/3 = 10, m>A is at least 2.5 times larger than the LSP mass, and could be arbitrarily large 
if very large values of m are included in the scan. On the other hand, for tan /? = 50, only 
a very narrow range of m^ values is possible for a given LSP mass, indicating that now 
has become almost (although not quite) independent of m . 

The (m^o,rrih) plane is shown in Fig. 11. In the left frame we see that the upper black 
(DM-allowed) region is also almost horizontal, indicating that in this "focus point" region 
rrih depends very weakly on mi/2- In fact, since we have m\ ^> m\ everywhere in this frame, 
rrih depends on mo and m\/2 only through loop corrections, in particular through the soft 
breaking stop masses, which are not sensitive to mi/2 as long as m\ ^> m\i 2 - Since in Fig. 1 
this region covers only a narrow range of m$ it also corresponds to a narrow range of m^. On 
the other hand, in the f\ co-annihilation region we have m\ <C m\ / 2 , so that the t masses are 
mostly determined by mu2- Since this region covers a broad range of m-y/2, it also extends 
over a significant range of m^. In this second branch the logarithmic increase of m^ with 
the SUSY mass scale, here represented by the LSP mass, is clearly visible. 

For tan j3 = 50 the four distinct DM-allowed bands shown in the corresponding Figs. 2 
or 7 collapse to two bands. The fi co-annihilation band is connected to the band where 
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Figure 9: The mSUGRA (mi ,rrif 1 ) parameter space with all constraints imposed for Aq = 0, 
H > 0,tan/3 = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 
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Figure 10: The mSUGRA (m^o, uia) parameter space with all constraints imposed for A = 0, 
fx > 0,tan/3 = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 



2m^o is slightly above tjia even in the (mo, m^) plane. Together with the band where 2m^o 
is slightly below m^, which gives very similar results for m^, they form the lower black band 
in Fig. 11. The DM-allowed region with mixed higgsino-bino LSP, which is quite prominent 
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Figure 11: The mSUGRA {m^o, rrih) parameter space with all constraints imposed for Aq = 0, 
fi > 0,tan/? = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 

for tan/5 = 50, gives the upper black band in the right frame of Fig. 11. Note also that for 
fixed m^o, in the DM-allowed region for tan/? = 50 exceeds that for tan/? = 10 only by 
about 1 GeV. However, in the latter case m^o up to ~ 700 GeV can be compatible with the 
DM constraint, whereas for tan/? = 10 this requires m^o < 500 GeV; the upper bound on 
rrih in the DM-allowed region therefore grows by about three GeV when tan/3 is increased 
from 10 to 50. These results are compatible with those of ref. |20j . 

Finally, the (mh,mA) plane is shown in Fig. 12. The most obvious feature here is the 
strong correlation of these two Higgs boson masses. In the left frame (tan/3 = 10) the lower 
black (DM-allowed) strip is the f\ co-annihilation region, whereas the upper black strip 
is the "focus point" region. In the latter part of parameter space is relatively small, 
and effects of the bottom Yukawa coupling are still almost negligible, so that tua — itlq. 
Moreover, both and are quite insensitive to mi/2, so long as 3> rriy 2 . This leads 
to a strong compression of the accessible region, which in this part of parameter space almost 
coincides with the DM-allowed region. However, the upper black strip in the right frame is 
surrounded by "inaccessible" light grey regions only because we limited our scan to m < 3 
TeV; otherwise it would be connected by the (experimentally and theoretically allowed, but 
DM-disfavored) white region. 

In the right frame of Fig. 12, i.e. for tan/? = 50, the lower black region is the "focus 
point" region, whereas the upper black strip is the f\ co-annihilation region merged with 
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Figure 12: The mSUGRA (rrih, vha) 'parameter space with all constraints imposed for Aq = 0, 
H > 0, tan/3 = 10 (left) and tan/3 = 50 (right). The top mass is fixed to m t = 172.7 GeV. 

the A— pole region, as discussed in the context of Fig. 11. For given itla, rn^ is maximized 
if mi/2 is minimized, i.e. in the "focus point" region. The reason is that rriA is sensitive to 
mi/2 already at the tree-level, through the relation that fixes /i 2 in terms of M| and the 
soft breaking parameters; in contrast, if m\ Mf , rrih depends on 7771/2 only through loop 
effects. We also note that the accessible region of the (777^,777^4) plane becomes very narrow 
for large tan (3. 

4. Sparticle and Higgs boson mass bounds 

The figures shown in Sec. 3.1 show that the allowed region in the (7770,7711/2) plane depends 
very strongly on the value of tan/3. There is also a significant dependence on A . Finally, 
even though the top mass is now the (relatively) most accurately known of all quark masses, 
we saw that varying m t within its current limits still moves the boundaries of allowed regions 
by hundreds or, in case of the "focus point" region, even thousands of GeV. Similar shifts 
of the allowed region occur when plotted in terms of pairs of physical masses, as shown in 
Sec. 3.2. 

One simple way to show the total allowed ranges of physical masses is to scan over the 
entire parameter space that is consistent with a given set of constraints; this is the topic of 
this Section. We saw in Sec. 2 that not all constraints should be treated on an equal footing. 
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Briefly, lower bounds on masses (or cross sections or branching ratios) from accelerator-based 
experiments are most robust, since both beam and detector are (hopefully) well controlled 
by the experimenter. Bounds on masses, and on cross sections of processes that can occur 
at tree level, usually also do not have many theoretical ambiguities. In contrast, we saw 
that one can evade the b — > 57 constraint by a relatively minor modification of the model 
|TS]. In the case of the <? M — 2, there is the additional ambiguity due to the ~ 2a discrepancy 
between SM 'predictions' based on different data sets, see eqs. (dJ) and ©• Finally, the 
DM constraint (J7|) required several (reasonable) assumptions for its derivation, and needs 
additional assumptions to be translated into allowed regions in the mSUGRA parameter 
space. 

It was originally hoped that the upper bound on the DM relic density (the so-called 
'overclosure' constraint) would allow to establish reliable, useful upper bounds on sparticle 
masses. Under the standard assumptions listed in Sec. 2, the predicted LSP relic density is 
proportional to the inverse of the (effective) LSP annihilation cross section into SM particles 
(or MSSM Higgs bosons, if kinematically accessible). This cross section in turn (through 
dimensional arguments, or by unitarity) scales like the inverse square of the relevant mass 
scale. Indeed, unitarity does allow to establish an upper bound on the mass of any WIMP; 
however, this bound exceeds 100 TeV |5()j . and is thus not particularly useful, since we lack 
the means to build colliders that could cover this kind of mass range. In the context of 
mSUGRA, it became clear quite early on that very, even "unnaturally", large masses 
can be compatible with the DM constraint (jZJ even in standard cosmology. 

On the other hand, we did see in Figs. 1-6 that this constraint excludes large chunks of 
otherwise allowed parameter space. One might therefore think that it would at least affect 
the lower bounds on sparticle masses significantly. Table 1 shows that this is not really 
the case. This table lists lower bounds on the masses of some new (s)particles, imposing 
various sets of constraints. We always impose all constraints discussed up to and including 
eq. (JIJ in Sec. 2. Sets I through III and IV through VI differ in that they are based on the 
conservative g^ — 2 constraint (jlj) and the more aggressive constraint (J5jl . respectively. Since 
the latter requires a positive supersymmetric contribution to g^, it allows us to derive upper 
bounds as well as lower bounds on the masses of new (s)particles. In addition to these basic 
constraints, sets II and V impose the b — > 57 constraint (JHJ), and sets III and VI in addition 
impose the DM constraint (jZJ. 

We see that the lower bounds on the masses of some key (s)particles always saturate 
their current bounds from collider physics, no matter what additional constraints we impose. 
This is true, in particular, for the lighter chargino, the lightest charged slepton (always f\ 
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Table 1: Lower bounds on the masses of superparticles and Higgs bosons, and upper bound 
on the LSP-nucleon scattering cross section, derived in mSUGRA under six different sets 
of assumptions. "HWIP" and "HSIP" stand for "heaviest weakly interacting particle" and 
"heaviest strongly interacting particle", respectively. In all cases the constraints discussed in 
Sec. 2 up to the requirement are imposed. In Set II we in addition impose the constraint 
(QJ) from b — > 57 decays (including the sign of the decay amplitude). Including in addition 
the DM constraint |7p leads to Set III. Sets IV-VI are like Sets I-III, except that the more 
conservative g^ — 1 constraint has been replaced by the more aggressive requirement 
in these cases we give allowed ranges, rather than only one-sided bounds. (Upper bounds on 
masses in Set V are the same as in Set IV.) All limits have been obtained by scanning the 
parameter space (GP, for 166.9 GeV < m t < 178.5 GeV. 



Quantity 


Set I 


Set II 


Set III 


Set IV 


Set V 


Set VI 


rrie R ^ m- m [GeV] 


106 


106 


107 


[106,1320] 


106 


[108,1300] 


m~ eL ~ m AL [GeV] 


152 


168 


169 


[152,1330] 


168 


[171,1310] 


m fl [GeV] 


99 


99 


99 


[99,1020] 


99 


[99,915] 


rrif 2 [GeV] 


156 


171 


172 


[156,1160] 


171 


[174,1130] 


m, T [GeV] 


130 


149 


149 


[130,1160] 


149 


[152,1120] 


[GeV] 


105 


105 


105 


[105,674] 


105 


[105,667] 


m-± [GeV] 


218 


218 


233 


[219,1003] 


227 


[337,999] 


m^o [GeV] 


52 


52 


53 


[52,359] 


53 


[55,357] 


m^o [GeV] 


105 


105 


105 


[105,674] 


105 


[106,667] 


m^o [GeV] 


135 


135 


135 


[135,996] 


135 


[292,991] 


m^o [GeV] 


217 


218 


234 


[218,1003] 


226 


[337,999] 


m- g [GeV] 


359 


380 


380 


[361,1880] 


399 


[412,1870] 


m d R - m s R [GeV] 


406 


498 


498 


[406,1740] 


498 


[498,1740] 


m d L - m 'sL [GeV] 


424 


518 


518 


[424,1810] 


518 


[518,1800] 


m- bi [GeV] 


294 


459 


463 


[295,1520] 


459 


[463,1500] 


m~ b2 [GeV] 


400 


498 


498 


[400,1600] 


498 


[498,1590] 


m k [GeV] 


102 


104 


104 


[102,1440] 


231 


[244,1440] 


m h [GeV] 


429 


547 


547 


[431,1600] 


547 


[547,1590] 


m h [GeV] 


91 


91 


91 


[91,124] 


91 


[91,124] 


m H [GeV] 


111 


111 


111 


[111,975] 


111 


[111,954] 


m H ± [GeV] 


128 


128 


128 


[128,979] 


128 


[128,960] 


rnnwip [GeV] 


349 


362 


366 


[351,1330] 


366 


[371,1310] 


™hsip [GeV] 


432 


556 


556 


[432,1880] 


556 


[566,1870] 


(JyOp [ab] 


140 


140 


7.5 


[10- 4 ,140] 


140 


[10- 4 ,7.5] 



in mSUGRA), as well as both the CP-even Higgs bosons. The lower bound on the masses 
of the charged Higgs boson is also independent of the constraints imposed; it follows almost 
directly from the structure of the MSSM. The lower bounds on the masses of the gluino and 
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the lighter two neutralinos are to a large extent fixed by gaugino mass unification and the 
chargino mass bound. In particular, the bound on m^o is practically identical to that on 
m^± . Note that gaugino mass unification holds for running (DR) masses. Going to the pole 
mass entails potentially quite substantial radiative corrections in case of the gluino [51]. The 
lower bound on rrig therefore increases by about 15% when going from the loosest constraints 
(Set I) to the tightest ones (Set VI), the biggest increase being due to the b — > 37 constraint, 
which excludes scenarios with large tan (3 and relatively small squark masses. 

The same effect is also visible in the lower bounds on first and second generation squark 
masses themselves, which increase by about 20% when this constraint is imposed. Note 
that even for the loosest set of constraints, Set I, the lower bounds on the masses of first 
and second generation squarks are significantly above current search limits [2Z|- This is a 
consequence of the assumed universality of scalar masses, together with the requirement of 
having sufficiently large soft breaking masses in the stop sector to satisfy the Higgs search 
limits. These bounds are therefore saturated for the largest possible m t value. 

The masses of the heavier neutralinos and charginos also depend only weakly on the set 
of constraints imposed. Most of these lower bounds again follow directly from the structure 
of the MSSM with gaugino mass unification, together with the LEP bound on m~±; the 
universality of scalar masses and A— parameters, which are defining properties of mSUGRA, 
play little role here. The only significant exception is the increase of m^o for the most 
restrictive set of constraints (Set VI). This bound is saturated in the "focus point" region of 
large mo, where the supersymmetric contribution to — 2 tends to be below the range (|5f . 

The lower bounds on the masses of third generation squarks are the quantities that are 
most sensitive to the additional constraints imposed in Sets II through VI. In particular, 
requiring both the more aggressive — 2 constraint (j5]) and the b — > 57 constraint (jUJ) 
more than doubles the lower bound on m$ 1 , from the LEP limit of ~ 100 GeV that can 
be saturated for Sets I through IV, to about 240 GeV. Combinations of parameters leading 
to a light ti which are compatible with the b —>■ sj constraint have relatively small tan/?, 
but large mo; this leads to a supersymmetric contribution to — 2 below the range ([5]). 
One can also generate a light t\ by taking modest values of mo and m.1/2, in agreement with 
the aggressive g^ — 1 constraint (J5J); however, this requires very large values of |Ao|/mo, 
which leads to a violation of the b — > S7 constraint (JHJ). This latter constraint by itself also 
increases the lower bound on m~ bi by about 50%, since a light b\ requires large tan/3 (which 
maximizes the bottom Yukawa coupling, as well as bi — &r mixing), which in turn leads to 
large (negative, for fi > 0) supersymmetric contributions to b — > 57 decays. 

In contrast, imposing in addition the DM constraint (J7J) has very little effect on the lower 



23 



bounds on sparticle masses. It does, however, drastically reduce the maximal possible elastic 
spin-independent LSP-proton scattering cross section, which is shown in the last line of the 
Table. The calculation of this cross section is based on refs. [22] • I n mSUGRA the potentially 
largest contributions come from the exchange of the heavier CP-even Higgs boson. Since 
increasing tan/? both reduces the mass of this boson and increases its coupling to down- 
type quarks, the cross section grows quickly when tan/3 becomes larger. In addition, it is 
maximized by increasing the coupling of the LSP to Higgs bosons, which requires rather 
larger gaugino-higgsino mixing; this cross section is therefore largest in the "focus point" 
region However, the same coupling also leads to too effective LSP annihilation, resulting 
in too low a relic density. Imposing the lower bound on the relic density in (JJJ) therefore 
reduces the upper bound on this cross section by about a factor of 20. We should mention 
here that even the reduced value of 7.5 ab, which is saturated for an LSP mass near 160 GeV, 
exceeds the current experimental lower limit on this quantity, if standard assumptions about 
the distribution of the LSPs in the halo of our galaxy are correct jHl] • Due to the uncertainties 
in this distribution we have not included LSP search limits in our set of constraints. 

While all sets of constraints allow the masses of some weakly interacting sparticles to lie 
right at the current experimental limit, mSUGRA implies that the heaviest weakly interacting 
new particle (sparticle or Higgs boson) must lie above ~ 350 GeV at least. Note that this 
limit lies well above the lower bound on the mass of any one weakly interacting (s)particle. 
The reason is that these bounds cannot be saturated simultaneously. For example, the lower 
bounds on slepton masses are saturated at moderate values of tan/3, big enough to avoid 
excessive lower bounds from Higgs searches, but not so big as to imply strong lower bounds 
from the — 2 constraint (@}. In contrast, the lower bounds on the masses of the heavy 
Higgs bosons are saturated at very large values of tan/3. Similarly, the lower bound on the 
mass of the heaviest strongly interacting sparticle in a given spectrum is somewhat larger 
than the largest lower bound considered separately. 

As mentioned earlier, imposing the aggressive — 2 constraint (J5J) allows to derive upper 
bounds in the masses of sparticles and Higgs bosons. The reason is that the supersymmetric 
contribution, which comes from gaugino-slepton loops, would vanish if either the gaugino 
masses or the slepton masses became very large. This leads to upper bounds on both mo 
and mi/2, as can be seen by studying the blue regions in Figs. 1 through 6. This in turn 
imposes an upper bound on via the condition of radiative electroweak symmetry breaking. 
As a result, all sparticle and Higgs masses can be bounded in mSUGRA using this single 
constraint! We emphasize that one needs to assume universality of both scalar and gaugino 
masses to derive these constraints. Numerically, the upper bounds on the masses of the first 
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generation squarks as well as on the gluino mass imply that discovery of supersymmetry 
at the LHC should be straightforward [22] ■ Unfortunately even in this favorable scenario 
discovery of charginos would not be guaranteed even at a 1 TeV e + e~ collider, and discovery 
of sleptons would be guaranteed only at a CLIC-like machine operating at s/s > 3 TeV. 
On the other hand, imposing this constraint reduces the upper bound on to 124 GeV 
(which becomes 127 GeV if one allows for a 3 GeV theoretical uncertainty), which might 
perhaps even give the Tevatron a chance to probe this scenario (with probably less than 
compelling statistical significance, alas) [HE] • Imposing in addition the b — > 57 constraint 
does not change these upper bounds at all; even imposing the DM constraint (J7J) leaves these 
upper bounds almost unchanged. 

One lesson from Table 1 is that imposing the DM constraint (J2J) has little effect on either 
upper or lower bounds on the masses of sparticles and Higgs bosons in mSUGRA, if one 
scans over the entire allowed parameter space. Figs. 13 show that the dramatic reduction of 
the allowed region that results from this constraint that is evident in the Figures presented 
in Sec. 3 does narrow down the possible ranges of masses when tan/3 is kept fixed. These 
figures compare the allowed ranges of the masses of e#, xf , h an d H for the constraint Sets 
IV and VI of Table 1. In the former case one can saturate the LEP limit on the mass of 
the lighter chargino for any tan (3 > 5. However, light charginos are DM-compatible only 
in the "focus point" region, which in turn is (barely) compatible with the — 2 constraint 

only at large tan j3. Moreover, at very large tan /3 the b — > 57 constraint becomes quite 
severe. As a result, constraint Set VI allows to saturate the LEP chargino mass bound only 
for 40 < tan/3 < 50. 

The combined effect of the DM and b — > 37 constraints on the lower bound on mi is even 
more dramatic. Without these constraints, the LEP limit on this mass can be saturated for 
any tan/3 < 50. However, as we already saw in Table 1, the b — > 57 constraint increases the 
lower bound on this mass bound by more than a factor of 2, if one insists on the aggressive 
Qfj, — 2 constraint (J2J); Fig. 13 shows that this lower bound then also increases quite rapidly 
with increasing tan/3. As a result, if tan/3 was known, imposing constraint Set VI would 
allow to predict to within a factor of ~ 3. However, since the allowed band moves 
upward with increasing tan/3, we can currently predict only within a factor of ~ 6, even 
if we impose this most restrictive set of constraints, as shown in Table 1. 

Imposing the DM constraint (|ZJ) reduces the upper bounds on sparticle masses for fixed 
tan/3. This effect is quite mild in most cases, but becomes prominent for cr at small and 
moderate values of tan/3. The g M — 2 constraint (J2J) by itself already leads to a strong tan/3 
dependence of this bound; recall that the supersymmetric contribution to this quantity is 
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Set IV: g -2 from e e data, no DM constraint 
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Set VI: g -2 from e e data, including b -> sj and DM constraints 
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Figure 13: The minimal and maximal value of select sparticle and Higgs masses: solid (black) 
curves: e^; dashed (red) curves: t\; dotted (green) curves: xf; dot-dashed (blue) curves: H. 
The upper (lower) figure shows the bounds without (with) imposing the DM and b — > 57 
constraints. 
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essentially proportional to tan/5. We saw at the end of Sec. 3.1 that for small and moderate 
tan/? the only overlap of the DM and (aggressive) — 2 allowed regions occurs in the f 
co-annihilation region, which has relatively large gaugino masses; as a result, one needs even 
smaller slepton masses to produce a sufficiently large g M — 2. However, once tan j3 > 40, one 
can satisfy this last constraint even in the "focus point" region; imposing the DM constraint 
then has little effect on the upper bound on rrie R - Finally, we note that the allowed range 
of rriii is fixed almost completely by the "base" set of constraints plus the aggressive — 2 
constraint (jSJ); imposing in addition the b — > sj and DM constraints has little effect here. 

5. Summary and Conclusions 

In this paper we provide an updated scan of the mSUGRA parameter space. This includes 
the new central value of the mass of the top quark, the inclusion of additional higher order 
corrections to the mass of the lightest CP-even Higgs boson h, and new information on the 
sign of the matrix element for b — > 57 decays from the analysis of b — > s£ + £~ decays. 

We find that the reduction of the central value of m t from 178 GeV to about 173 GeV 
shifts the allowed parameter space significantly. This is the consequence of two effects: 

• The corrections to m\ scale like the fourth power of m t , but only scale logarithmically 
with the sparticle (mostly stop) mass scale. As a result, a few percent reduction of 
m t has to be compensated by an increase of mi of up to several tens of percent. This 
relative shift increases with mj, and is therefore most prominent for smaller tan/3, 
where the LEP Higgs search constraint is most severe. 

• The location of the "focus point" region at m\ ^> m\ r 2 where the LSP acquires a sig- 
nificant higgsino component depends very sensitively on m t . The higgsino component 
of Xi is only sizable if is rather small. This will happen if the squared soft breaking 
mass of the Higgs boson coupling to the top quark is small or positive at the "weak" 
scale Qw, which should be chosen close to -/m^m^ to get good convergence of the 
perturbative series. This parameter in turn can be written as 

m H u (Qw) = am l + fom i/2 + cA l + dA m 1/2 , (8) 

where the dimensionless coefficients a, b, c, d depend on the dimensionless MSSM cou- 
plings as well as (logarithmically) on Qw The crucial observation is that \a\ <C 1 for 
m t ~ 175 GeV and not too small tan/3 (recall that the top Yukawa coupling at the 
weak scale is oc 1/ sin/3). In contrast, |6| is quite sizable, with b < 0. A small tti 2 h u {Qw) 
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is therefore only possible if ml ^> m\, 2 and a > 0. Since a depends (roughly) quadrat- 
ically on m t / sin/5, but only logarithmically on the sparticle mass scale (through Qw), 
a small change of m t therefore leads to a large shift of the value of Qy/i or '> equivalently, 
of mo where the "focus point" region starts. 

The first effect makes it more difficult to reconcile, at low values of tan /3, the Higgs mass 
constraint from LEP with the evidence for a positive contribution to the anomalous magnetic 
moment of the muon. However, for larger tan f3 the Higgs mass bound allows smaller sparticle 
masses, while the contribution to — 2 remains significant for larger sparticle masses. As a 
result, even for m t = 173 GeV both constraints can be satisfied simultaneously for tan (3 > 10. 

By far the most stringent constraint on the parameter space comes from the requirement 
that the lightest neutralino should have the correct relic density. As discussed in Sec. 2, this 
constraint can only be translated into bounds on the mSUGRA parameter space if several 
assumptions are made. Under the usual assumption of thermal LSP production and standard 
cosmology, only a few discrete "DM allowed" regions survive. Out of these, the "bulk" region 
of moderate mo and moderate rai/2 is affected most by the reduced mass of the top quark; in 
fact, it disappears altogether if m t is indeed near its current central value of ~ 173 GeV. The 
T\ co-annihilation region is also reduced in size, since the region excluded by Higgs searches 
now extends to larger values of my 2- On the other hand, as noted above, the region where 
the LSP has significant higgsino component becomes larger when m t is reduced. Similarly, 
the lowest value of tan f3 where 2m^o ~ m& (the A— pole region) becomes smaller; this region 
also becomes broader, thanks to the increased strength of the AxiXi coupling. In contrast, 
the conceptually similar h— pole region is much reduced in size. 

We also provided views of the mSUGRA parameter space plotted in the plane spanned 
by two physical sparticle or Higgs boson masses. Whereas some masses are essentially 
independent of each other (e.g. mg R and m%o), others are strongly correlated (e.g. and 
nth); most pairs are intermediate between these extremes. 

In addition to these scans of parameter space, we also provided upper and lower bounds 
on the masses of Higgs bosons and sparticles. Here it is crucial to properly include the 
uncertainty of the input parameters, in particular, of m t . This sensitivity comes mostly 
through the dimensionless coefficients in eq.(jSJ), as well as the analogous expression for m\ . 
These coefficients determine through the condition of electroweak symmetry breaking, 
which affects the spectra of neutralinos, charginos and Higgs bosons. The latter is also 
directly dependent on these coefficients; through the Higgs search limits, they then affect 
(the lower bounds on) all sparticle masses. It should be noted that these coefficients also 
depend on other input parameters, in particular on as and m^. However, this dependence is 
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much milder than that on m t . We therefore believe that varying m t over its entire currently 
allowed 2a range, while keeping a 5 and m h fixed to their default values, gives a reasonable 
estimate of the effect of the current input parameter uncertainties. 

The somewhat surprising result of these scans is that the masses of many superparti- 
cles and Higgs bosons can still lie right at their current limits from collider searches even 
if the most restrictive set of constraints is applied, including the Dark Matter constraint 
and the more aggressive interpretation of the — 2 constraint. This means that ongoing 
and near-future experiments still have good chances to discover new particles even in this 
very constrained version of the MSSM. Not all these lower bounds can be saturated simul- 
taneously, however. As a result, the most robust constraints (essentially the collider limits 
plus a conservative version of the g^ — 2 constraint, with no DM requirement) by themselves 
already imply that a 500 GeV linear collider will not be able to discover all new weakly in- 
teracting (s)particles; one will need an energy of at least ~ 700 GeV to achieve this. We also 
saw that these lower bounds are in most cases surprisingly insensitive to the introduction 
of new constraints; in particular, requiring the lightest neutralino to be a good thermal DM 
candidate does not shift them much. 

Useful upper bounds on the masses of sparticles and Higgs bosons (with the exception 
of the lightest CP-even Higgs boson, see ref. can only be derived if we assume that 
a positive supersymmetric contribution to g M — 2 is indeed required, as is indicated (at the 
~ 2.5a level) when data from e + e~ — > hadrons are used to calculate the SM prediction for 
this quantity. This imposes upper bounds on the masses of both sleptons and gauginos; 
in the mSUGRA context this implies upper bounds on both mo and mi/2, which leads to 
upper bounds on all new (s)particles. Quantitatively, we find that this requirement by itself 
implies that strongly interacting sparticles must be within the reach of the LHC. Moreover, 
a ~ 1 TeV e + e~ collider would then be able to at least discover superparticles, in the X1X2 
channel. However, even in this case one may need a CLIC-like collider, with center of mass 
energy nearly reaching 3 TeV, to discover all new weakly interacting (s)particles. We stress 
again that these (upper) bounds have been obtained by scanning over the entire parameter 
space, including scanning over the 2a range of m t . Imposing in addition the Dark Matter 
constraint narrows down the allowed ranges of some masses if tan/5 is held fixed, but has 
little effect on the absolute upper bounds after scanning over the entire allowed parameter 
space. 

We conclude that mSUGRA remains viable. In fact, even after imposing all plausible 
experimental and theoretical constraints the allowed parameter space still contains a large 
variety of different spectra, with quite different phenomenology. Even if supersymmetry 
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provides the missing Dark Matter in the Universe and explains the possible excess of the 
anomalous dipole moment of the muon, superparticles might be out of reach of both the 
Tevatron and the first stage of a linear e + e~ collider; on the other hand, it is also entirely 
possible that their discovery is "just around the corner" . 
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